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Abstract 

^ I A self-consistent Quasiparticle-Random-Phase- Approximation (QRPA) model which employs the 

canonical Hartree-Fock-Bogoliubov (HFB) basis and an energy-density functional with a Skyrme 
mean field part and a density-dependent pairing, is used to study the monopole collective excita- 
Q^ ', tions of spherical even-even nuclei. The influence of the spurious state on the strength function 

o 

I of the isoscalar monopole excitations is clearly assessed. We compare the effect of different kinds 

of pairing forces (volume pairing, surface pairing and mixed pairing) on the monopole excita- 

OO . 

. tion strength function. The energy of the Isoscalar Giant Monopole Resonance (ISGMR), which 

is related to the nuclear incompressibility Koo, is calculated for tin isotopes and the results are 
^ ■ discussed. 
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I. INTRODUCTION 



The nuclear structure community have made great achievements in understanding the 
structure of the ground state and of the excited states of stable atomic nuclei. As radioactive 
beams provide more experimental results on the nuclei far from the stability valley, the 
challenge is how to extrapolate the theoretical models and predict or describe in detail the 
exotic properties of the nuclei with large neutron or proton excess. Another challenge is the 
prediction of the properties of nuclear matter in a broad range of densities, i.e., in connection 
with neutron stars, and to understand the origin of these new properties. 

For medium-mass and heavy nuclei, the most microscopic models that we can use are 
the mean field models based on the effective interactions, either in the non-relativistic or 
relativistic framework. For closed shell nuclei, Hartree-Fock (HF) theory has already been 
proven to be a powerful tool to describe thej)roperties of ground states in particular 



using the zero-range Skyrme interactions 



In the open shell nuclei, the effect 



of nuclear pairing shows up. A simple theory for the ground-state pairing is HF+Bardeen- 



Cooper-Schriffer (BCS)[3]. The nuclei close to the neutron or proton drip lines, may exhibit 
some very unusual features such as the neutron or proton skin 8| and the neutron haloes ^. 
In these very neutron-rich or proton-rich nuclei, nuclear pairing plays a crucial role for the 
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14l |. A more appropriate 



theoretical understanding of these new phenomena 

n 

theory is the HFB approach [1[ because the pairing component can no longer be treated 
as a residual interaction, i.e., a small perturbation important only in the neighborhood of 
the Fermi surface, as in the nuclei close to the line of beta stability. This is seen from 
the approximate HFB relations between the Fermi level A, pairing gap A, and the particle 
separation energy s: s ~ —A — A, since s is very small and A + A ^ 0, for drip-line 
nuclei. Consequently, the mean field characterized by A and the pairing field A are equally 
important. Using appropriate effective interactions in mean field and pairing field, the 
HFB approach is already sophisticated enough to allow precise analysis of ground state 
properties, i.e., binding energies, average neutron pairing gaps, etc., in most nuclei either 
using the Skyrme force plus a density-dependent pairing force |l5] or the finite-range pairing 



force 
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17|. 



Studying the nuclear collective excitations is another important tool to understand the 
structure of nuclei and predict the exotic properties of nuclei far from stability valley or the 
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properties of nuclear matter. The QRPA is a standard method for describing these collective 
excitations in open-shell superconducting nuclei with stable mean field solutions 
Important nuclear collective excitations are the nuclear compressional modes - particularly 
the ISGMR - which provide the optimal route to determine the nuclear incompressibility [l9|, 
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2l|. Both non-relativistic RPA 
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2^ and relativistic RPA or QRPA |25|, l26| were 



recently used in studying the nuclear collective excitations and the nuclear incompressibility. 
However, the fully self-consistent QRPA, formulated in the HFB canonical basis, which 
was introduced and accurately tested using Skyrme energy density functionals and density- 
dependent pairing functionals in Ref. 27|] has not been applied to extract information about 
the nuclear incompressibility. If the models are characterized by a nuclear incompressibility 
around 230~240 MeV (or 250~270 MeV in the case of relativistic mean field (RMF)), 
they will give the right ISGMR centroid energies compared with the experimental data in 



208 



Pb 
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25| . However, with these modes one tends to overestimate the centroid energies 



in the tin isotopes 
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261]. Is there any need of more sophisticated considerations for the 



nuclear pairing in tin ? To answer this question, we study the ground state properties of 
tin isotopes using the HFB theory with a Skyrme mean field and different kind of pairing 
forces, and then we apply a self-consistent QRPA model for the ISGMR. 

The paper is organized as follows. Section II discusses the QRPA employed in this study. 
Since the QRPA equation is solved on the canonical HFB basis, in this section we also 
present a survey of the canonical HFB theory. We mention briefly the meaning and role of 
self-consistency in the context of RPA and QRPA calculations. In section HI, we give the 
numerical details and the results. In this section, we discuss the method of projecting out the 
spurious state, the influence of pairing force on the strength function, and in particular, we 
study the energy of the ISGMR for tin isotopes. The last section contains a short summary 
and our conclusions. 



II. QRPA WITH CANONICAL HFB BASIS 

The standard approach for describing collective excitations in magic nuclei is the so-called 
Random-Phase-Approximation (RPA) with stable mean-field solutions. The extension to 
open-shell nuclei is the QRPA which includes not only the particle-hole channel (mean 
field), but also the particle-particle channel (pairing field). The first step in solving the 
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self-consistent QRPA is to obtain the ground state properties by solving the HFB equations. 
Since the HFB theory is discussed in details in many textbooks and articles Q, , here we 
just present briefly the HFB theory in canonical basis. 

The canonical basis is a particular basis in which the density matrix is diagonal, and the 
canonical basis wave functions are all localized. Starting from the effective Hamiltonian of 
a nuclear system, using the variational principle, the HFB equations in the canonical basis 



are obtained as 



15| 



{h-\)^u{u^v^ + UyV^,) + h^,y{u^,Uy-VyV^,) = (1) 
{h - \)t,y{Uf,Uy - VyVfj) - hfj,u{u^,Vy + UyV^,) = E^^, (2) 

where h and h are the particle-hole and particle-particle Hamiltonian in the canonical ba- 
sis, respectively. A is the Fermi energy while E^^u represents the HFB Hamiltonian in the 
canonical basis. The canonical states fi are not eigenstates of the HFB Hamiltonian. The 
right way to obtain the canonical basis is to solve the HFB equations in the quasiparti- 
cle basis to obtain the quasiparticle states, i.e., the quasiparticle wave functions (including 
the upper and lower conponents (pi and ^2) and quasiparticle ener gies , and then to build 
the density matrix and to diagonalize it to get the canonical basis 15]. The square of 
(m^) represent the probability that a certain state /i is occupied (empty). These quantities 
are solely determined by the diagonal matrix elements of particle-hole and particle-particle 
Hamiltonians, 




In the case of spherical symmetry, and without mixing of the proton and neutron states, 
the quasiparticle wave functions have the good quantum numbers {nljm); n turns out to 
be a good quantum number because the continuum is discretized inside a spherical box. 
Furthermore, the radial part of the quasiparticle wave functions can be chosen to be real. 
Then the quasiparticle wave function can be written as: 

Mnlj, ra) = "^Mllll ^ Y^]{r)^:^]{lmi\m,\jm) , z = 1, 2, (4) 

minis 

where i = 1,2 label the upper and lower components of the wave functions. The density 
matrix elements can be constructed using the lower components of the quasiparticle wave 
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functions, 



P(0-)(r,r') = ^02Nj,r)0;(n/j,r'), 



(5) 



where the sum runs over all the states for a given (/, j) block. We will get the canonical 
basis after diagonalizing the density matrix. 

The QRPA takes the quasiparticle vacuum as the approximate ground state and is aimed 
at the description of small amplitude, collective excitations. Defining the phonon excitation 
operators in the angular momentum coupled representation, 



a>l3 
a>f3 



where 



a/1 + 

1 



1 + 5, 



a/3 



{jamajpmp\JM)alal 
{jamajpmp\JM)apaa 



(6) 
(7) 

(8) 

(9) 
(10) 



{-iy+''A^p{j-M), 

the QRPA equation can be easily obtained from the linearization of Schrodinger equation, 

^ j Baf3,^S | / \ ^ ^ / X^^ \ ^^^^ 

In the angular momentum coupled representation and canonical HFB basis, the matrix 



elements Aap^-^i and Bajs^-ys have the form 



m 



27|, 



A 



B. 



1 



(12) 



a/9,75 



(13) 



with 



Gap-^& = ^ {jamajf3mp\JM){j^m^jsm5\J'M')V^l^^ (14) 
Fap-yS = {jamajf3mp\JM){j^m^jsms\J'M')V^^^^, (15) 

where V^^^^ represent the matrix elements of the particle-hole effective interaction, and 
^aJ^-yS represent the matrix elements of the particle-particle effective interaction. Actu- 
ally, in the canonical basis, the equations of HFB+QRPA are almost the same as those of 
HF-BCS+QRPA, except for the presence of off-diagonal terms related to the quasiparticle 
energies in the former case. 

We remind the reader that self-consistency means, in the RPA context, that the particle- 
hole residual interaction is obtained as the exact second derivative of the energy functional 
with respect to the density. Up to quite recently, most of the RPA calculations were not fully 
self-consistent, rather some parts of the residual interaction were neglected. Only full self- 
consistency allows conservation of the energy-weighted sum rule (EWSR) according to the 



Thouless theorem 



28l | . Moreover, in the case of monopole the inclusion of the whole residual 



interaction is crucial in order to evaluate accurately the energy of the giant resonance 23| . 



and to assess quantitatively the value of Koo jl9| 



In the QRPA framework, self-consistency should mean, in addition to what has been 
stated for RPA, that the particle-particle residual interaction is obtained as the second 
derivative of the energy functional with respect to the anomalous density or pairing density. 
This prescription is obeyed both in HF-BCS+QRPA and in HFB+QRPA. In the former 
case, however, the model spaces for HF-BCS (which includes, as a rule, only one major 
shell above the Fermi energy) and for QRPA (which is much larger), differ. In this sense, 
only HFB+QRPA is, strictly speaking, self-consistent. The Thouless theorem concerning 

so valid for the self-consistent QRPA based 



29|. 



the EWSR obtained within HF plus RPA is a 
on HFB solution as it is demonstrated in Ref. 

The fully self-consistent calculation must include the pairing-rearrangement terms in 
the matrix elements, namely the contribution to the particle-hole interaction coming from 



the density dependence of the pairing force 27|, ISOj. In the case of the volume pairing 



force, there are consequently no pairing-rearrangement terms at all, while in the case of 
the surface pairing force there exists a contribution from the pairing-rearrangement terms. 



We have carefully checked that this contribution has a negligible effect on the ISGMR. 
For example, the ISGMR centroid energies obtained from the QRPA calculations based on 
canonical HFB solution without the effects of the pairing-rearrangement using the SKM* 
force and the surface pairing force are 16.32, 16.13, and 16.0 MeV for ^^^Sn, ^^^Sn and ^^^Sn, 
respectively. When including the effects of the pairing- rearrangement, the values of the 
ISGMR centroid energies change to 16.34, 16.11, and 16.02 MeV for ^^^Sn, ^^'^Sn and ^^^Sn, 
respectively. Therefore, in the figures and in the discussion below, we quote results obtained 
without the pairing-rearrangement terms. 

III. RESULTS AND DISCUSSIONS 

For the ground state properties of spherical nuclei, the HFB equations are solved in 
coordinate space in a spherical box. In the particle-hole channel, we use a Skyrme force, 
i.e., SLy5 |J] or SKM* [6|, and in the particle-particle channel, we use a zero-range density- 
dependent pairing force. 



where po = 0.16 fm"'^; the values of rj are 0, 1, and 0.5 for volume, surface and mixed pairing 
forces, respectively; the value of Vq is fixed by fitting the experimental data of the mean 
neutron gap of ^^°Sn (A„ = 1.321 MeV). In the case of the Skyrme force SLy5, Vo= -170.92, 
-537.05 and -272.99 MeV fm^ for volume, surface and mixed pairing forces, respectively. In 
the case of the Skyrme force SKM*, vo= -142.01, -490.48 and -233.22 MeV fm^ for volume, 
surface and mixed pairing forces, respectively. There are some numerical parameters in the 
actual calculations: (i) We use a high quasiparticle-energy cutoff (200 MeV) and a maximum 
angular momentum j^ax = 15/2. (ii) The radius of the box is fixed at 20 fm with a small 
mesh (0.05 fm). 

The self-consistency of the HFB+QRPA calculations requires the use in QRPA of a 
residual force derived from the HFB fields. At the same time, for a real self-consistent 
calculation, all the quasiparticle states produced by the HFB calculation must be used to 
build the matrices A and B in the QRPA equation. In actual calculation, there are about 
500-600 quasiparticle states in the HFB calculation up to 200 MeV which is the HFB cutoff. 
Some states with very small values of occupation probability in canonical basis or some 




(16) 
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two quasiparticle excitations with very high energy give httle contribution to the QRPA 
spectrum, so we cut off the canonical-basis wave functions by excluding those with very 
small values of occupation probability and with very high values of single particle energy, 
and check the convergence of the QRPA solution. Fig. [T] shows the isoscalar 0^ strength 
function (averaged with Lorentzians having 1 MeV width) in ^^°Sn with the Skyrme force 
SLy5 and the volume pairing force defined above, for different values of cutoff in excitation 
energy and occupation probability. Comparing the curve with circles with the curve with 
stars in Fig. [1], we find that the levels with occupation probabilities smaller than 10~^ have 
little influence on the isoscalar 0"^ strength function. Similarly, from the curves with squares, 
up-triangle, and circles in Fig. [H we find E^ut = 150 ~ 200 MeV is a suitable excitation 
energy cutoff in calculations to make our results stable at the level of 300 keV for the the 
centroid energy. Therefore, in the next calculations, the value of the cutoff on occupation 
probability is fixed at = 10~^, and the excitation energy cutoff is selected as E^ut = 200 
MeV. 



A. The spurious state 

In the spherical QRPA solution, there exists the problem of a spurious state in the 
monopole channel due to the particle number symmetry broken by the HFB solution 
The spurious state should be orthogonal to all other physical states, and should appear 
at zero energy. However, in actual calculations the spurious state is at low (but not zero) 
energy for small numerical inaccuracies, and contributes to the total strength. Therefore, 
the spurious state must be projected out from the real physical states. 

Starting from the actual quasiparticle QRPA set of states \n), we construct a new set of 
normalized states \n') which are the real physical states, 

\n') = \n) - an\s) (17) 

spurious state. The X- and F-amplitudes of the spurious state should be 



where \s) is the 
proportional to 



{ai3\N\0) = -2jaUaVa, (18) 
where \af3) is a pair of canonical states, N is the particle number operator, |0) is the phonon 
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vacuum state, and j = y/2j + 1. a„ is obtained by the condition (n'|A^|0) = 0, 



Or 



(19) 



Fig. [2] shows the influence of the spurious state on the isoscalar 0"*" strength function in 
^^°Sn with the Skyrme force SLy5 and the volume pairing force. The spurious state affects 
the low-lying region as it is expected, while it does not impact the ISGMR at all. 

B. The isoscalar 0+ mode 

The ISGMR is an important collective excitation because its excitation energy is re- 
lated to the nuclear incompressibility Kqo- The excitation operator of the ISGMR is 
^monopoie — SiLi ''^i' ^hc macroscopic picture for a state excited by this operator is the 
so-called "breathing mode" . The strength function is defined as 

SiE) = Y.\(^\P^'onopole\0)mE - (20) 

n 

and the moments of the strength function are 

ruk = j E^S{E)dE. (21) 

The centroid energy can be defined as the ratio between the EWSR sum rule mi and the 
non-energy weighted (NEWSR) sum rule itlq: Eq = mi/rrio. 

The ISGMR has been so far mainly described within the HF+RPA and the HF- 
BCS+QRPA which have been defined above. In Fig. [3] we show the results for the IS- 
GMR centroid energies in the tin isotopes using HF+RPA (by employing the Skyrme forces 
SLy5 [4], SGII Q|, and SKM* jol). The centroid energies are evaluated in the energy in- 



terval between 10.5 and 20.5 MeV, and the experimental data come from Ref. [31 



Other 



experiments devoted to the ISGMR in tin isotopes have been reported in Refs. [3J]. Their 
results are also shown, although the centroid energies are evaluated in the energy interval 
between 10 and 35 MeV. In the present work, we have no clear way to judge the difference 
between the two experimental results. The trend along the isotope chain is similar for the 
three forces. The force SLy5 {Koo = 233.8 MeV) gives higher centroid energy than the 
others (which have K^o around 215 MeV). HF+RPA overestimates the centroid energies by 
about 1 MeV in all the measured tin isotopes, and there exists a difference by few hundreds 



of keV between different Skyrme forces. The force SLy5 reproduces well the monopole en- 
ergy in ^°^Pb. The question why it overestimates this energy in the tin isotopes has been 
raised, as it is recalled in the introduction. If we take pairing into account, by means of the 
HF-BCS+QRPA, the results are very similar to those from HF+RPA, as it is expected and 
shown in Fig. HJ which presents the ISGMR centroid energies in the tin isotopes by using 
SLy5 and the volume pairing force. 

One possible guess could be that theoretical models overestimate the centroid energies in 
the tin isotopes because a more appropriate method to deal with the pairing correlations than 
HF-BCS+QRPA, i.e., the HFB+QRPA, is needed. As shown in Fig. ij HFB+QRPA makes 
the theoretical results significantly closer to experiment, especially in ii2,ii4,ii6gj^^ 

However, 

this improvement is not enough. The difference between HF-BCS+QRPA (HF+RPA) and 
HFB+QRPA comes from the fact that in ii2,ii4,ii6gj^^ ^^le single particle energies and occu- 
pation number of d3/2 and hi 1/2 are different for HF-BCS (HF) and HFB, while they are 
similar in ^^°'^^^'^^^Sn. These numbers are shown in Table [H 

Already from Fig. [3l it is clear that SKM* performs better than SLy5 in the tin isotopes. 
Having understood that pairing is also lowering the monopole energies, as a next step we 
study the isoscalar 0"*" excitations using HFB+QRPA with SKM* in the mean field. At the 
same time, we would also like to discuss in detail the influence of different kinds of pairing 
forces on the ISGMR centroid energies. Fig. [5] shows the isoscalar O"*" strength function 
(averaged with Lorentzians having 1 MeV width) in i20Sn with the Skyrme force SKM* for 
different kinds of pairing force. The results with the volume and the mixed pairing force 
are almost the same, while the results with the surface pairing force are different from the 
others, especially in the low energy region. This difference comes from the fact that in the 
ground state HFB calculations, the energies and occupation numbers of the levels around 
the Fermi surface are very similar for volume and mixed pairing, while they are different for 
surface pairing. Fig. [6] shows the ISGMR centroid energies in the tin isotopes with SKM* in 
the mean field and volume and surface pairing forces. Using the SKM* plus surface pairing 
force, one can reproduce the experimental data obtained at RCNP from ^^^Sn to ^^°Sn, while 
SKM* plus the volume pairing force still overestimates the the centroid energies by about 
0.5 MeV. The discrepancy with the experimental findings from Texas A& M remains quite 
significant. 
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IV. CONCLUSIONS 



The fully self-consistent QRPA based on the canonical HFB basis has been introduced by 
the authors of Ref . [2?^ . In the present work, we focus on the same model, having developed 
independently the formalism and the computer codes. Our motivation is the description 
of the isoscalar monopole excitations. The main novel aspect of this paper is the attempt 
to understand the anomaly in the experimental ISGMR energies of the Sn isotopes, which 
cannot be reproduced using effective forces which do reproduce the experimental energy in, 



26|, that is. 



e.g., ^^'^Pb; in other words, our goal is to answer the question raised in Ref. 
"why is tin so soft ?" . 

The model uses the zero-range Skyrme force in the particle-hole channel and the density- 
dependent pairing force in the particle-particle channel. The spurious state, caused by the 
fact that the particle number symmetry is broken within HFB and restored in QRPA, is 
projected out. We investigate the isoscalar O"*" strength function and the ISGMR centroid 
energies in the tin isotopes. We have also found that the pairing-rearrangement terms 
and find they have a negligible effects (about 0.1%) on the ISGMR centroid energies in 
the tin isotopes. Compared with HF-^RPA and HF-BCS-QRPA, HFB+QRPA makes the 
theoretical results of ISGMR centroid energies in the tin isotopes significantly closer to 
experiment. Pairing improves the absolute values of energy compared to experiment because 
of the attractive character of the particle-particle residual interaction. This amounts to say 
that the "softness" of tin with respect to monopole excitations is to some extent related to 
pairing. 

Using an appropriate Skyrme force like SKM* and surface pairing force, we reproduce the 
experimental ISGMR centroid energies from ^^^Sn to ^^°Sn. The pairing energy in nuclear 
matter, and consequently the effect of pairing on the nuclear incompressibility i^oo, is very 
small. Our results imply therefore that, whereas ^ospb leads to values of Koo around 230~240 
MeV, this value is about 10% smaller when Sn experimental data are used. Our results do 
solve partly the puzzle caused by the RCNP experimental results in the tin isotopes. At 
the same time, the 10% variation of from different nuclei points to our still incomplete 
understanding of the details of the nuclear effective functionals, in particular of their density 



dependence. Also in Fig. 8 of Ref. 



25| . it is evident that RMF calculations fit very well 



the energies of the monopole in Pb and Sn with two values of Koo which differ by ~10%. 
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Further investigations on how to constrain the surface and asymmetry contributions to the 
final nuclear incompressibility are demanded. 
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TABLE I: The single particle energies (MeV) and occupation number of d3/2 and hll/2 in the tin 
isotopes using the Skyrme force SLy5. In HFB, we use the volume pairing force in particle-particle 
channel. 





d3/2 


hll/2 




HF 




BCS 




HFB 




HF 




BCS 




HFB 




A 


e 




e 




e 




e 




e 




e 




112 


-8.69 





-8.63 


0.18 


-8.33 


0.21 


-7.14 





-6.90 


0.066 


-6.60 


0.048 


114 


-8.74 





-8.65 


0.26 


-8.37 


0.29 


-7.45 





-7.07 


0.089 


-6.80 


0.067 


116 


-8.69 





-8.67 


0.36 


-8.42 


0.40 


-7.49 





-7.23 


0.12 


-6.98 


0.095 


118 


-8.62 


0.5 


-8.70 


0.48 


-8.47 


0.54 


-7.61 





-7.35 


0.16 


-7.11 


0.14 


120 


-8.57 


1.0 


-8.74 


0.64 


-8.52 


0.68 


-7.72 





-7.46 


0.21 


-7.22 


0.20 


122 


-8.70 


1.0 


-8.81 


0.78 


-8.61 


0.80 


-7.69 


0.17 


-7.53 


0.29 


-7.28 


0.29 


124 


-8.83 


1.0 


-8.90 


0.87 


-8.71 


0.89 


-7.66 


0.33 


-7.57 


0.40 


-7.33 


0.41 
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FIG. 1: (Color online)The isoscalar O"*" strength function (averaged with Lorentzians having 1 MeV 
width) in ^^'^Sn using the Skyrme force SLy5 and the volume pairing force for different values of 
single particle energy cutoff and occupation probabilities. 
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FIG. 2: (Color online)The isoscalar O"*" strength function (averaged with Lorentzians having 1 
MeV width) including spurious state (line with circles) or not (line with squares) in ^■^''Sn using 
the Skyrme force SLy5 and the volume pairing force. 
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FIG. 3: Systematics of the moment ratios mi/rriQ for the ISGMR strength distributions in the tin 
isotopes using Skyrme HF+RPA. They are evaluated in the energy interval between 10.5 and 20.5 



MeV. The experimental data are extracted from Ref. 
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FIG. 4: Systematics of the moment ratios mi /mo for the ISGMR strength distributions in the tin 
isotopes. They are evaluated in the energy interval between 10.5 and 20.5 MeV. 
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FIG. 5: (Color online)The isoscalar O"*" strength function (averaged with Lorentzians having 1 MeV 
width) in i20sn with the Skyrme force SKM* for different kinds of pairing forces. 
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FIG. 6: Systematics of the moment ratios mi /mo for the ISGMR strength distributions in the tin 
isotopes using HFB+QRPA with the Skyrme force SKM* and volume and surface pairing forces. 
They are evaluated in the energy interval between 10.5 and 20.5 MeV. 
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